Comment on "Optimal probe wave function of weak-value amplification 
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First, a misconception about the spectrum of a confined particle is evidentiated. Then, the results 
are shown to be incorrect by means of a counter-example, an explicit preparation for the probe is 
given that yields an arbitrary amplification, and the source of the error is pointed out. 
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The authors of Rcf. [1] claim to have determined the 
optimal probe wave-function that provides the maxi- 
mum average readout for fixed preparation and post- 
selection. The results are independent of the actual cou- 
pling strength and they apply to operators satisfying 
A 2 = 1. For clarity and brevity we shall prepend the 
acronym SSH when referring to equations in Ref. [lj. 

After Eq. (SSH18), which we should take momentar- 
ily as good, Ref. [l[ states that "the position q of the 
probe takes the discrete value q — 2gn, with n being 
an integer because of the boundary condition in the mo- 
mentum space." This misleading statement is repeated 
in the caption of Fig. 2, where the function £i(q) is plot- 
ted at discrete points. This is the complementary form 
of the common misconception that the momentum of a 
confined particle is discrete. One should not confuse the 
indices of the Fourier basis for the expansion in energy 
eigenstates with the allowed values of the variable q or 
p. The position q has a continuous spectrum, even if 
the momentum is bounded to the interval [—ir/2g, w/2g\. 
To make this point more clear, let us consider a parti- 
cle of mass m confined in a one dimensional box with 
a coordinate 9 € [0, L]. Thus, the roles of q and p are 
temporarily switched in order to provide a familiar case. 
It is well known that the energy is quantized according 
to E n — p 2 l /2m, with p n — hnnj ' L. By analogy with the 
classical expression, it is tempting to identify p n with the 
allowed values of the particle's momentum. However, a 
von Neumann measurement of the operator P does not 
necessarily yield a value p n , but has a finite probability 
for any p. Even if the wavefunction is ip n oc sin {p n q), its 
Fourier transform ip n (p) gives a continuous probability 
distribution | ■(/>«. (p)| 2 that has two main peaks at p n and 
— p n but also several secondary peaks. 

Going back to the case studied in Ref. [lj, it is more 
convenient to work with the rescaled quantities x = q/g 
and k = gp, so that the pointer shift is (x) f — (q)f/g — 
(1 + |A u) | 2 )/2ReA ll ,. Then the Fourier transform of 
Eq. (SSH15) yields &(x) = i Q (x - (%), where 
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with a the sign of Re(A w ) 7 z = (1 — aA w )/(l+aA w ), and 
s, x) — J2k=o z I ( x + k) s the Lerch transcendent. In 



particular, at x = 2n, 
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for any integer n, while £0 (— a(2n + 1)) = 

(-z)V2|ReA,|/(l + crA w ), and £ (o-(2n + l)) = 
for n a positive integer. A further confirmation that 
£i should not be taken at discrete values is that 
E„ M 2n )\ 2 = 1/2, while Jdx\ii(x)\ 2 = 1, so that 
|£i(2n)| 2 do not represent probabilities. 

On the other hand, the claim that A(x) = (1 + 
|A tu | 2 )/2ReA MJ is the maximum pointer shift is un- 
founded. As a counter-example, let us consider a weak 
measurement of a spin 1/2 along the X-direction (no re- 
lation to the x representing the pointer readout), where 
the spin is initially prepared along the positive Z axis and 
afterward post-selected in a direction along the XZ plane 
forming an angle with Z . For an ideally weak measure- 
ment, obtained in the limit W — > of Eq. (SSH14), we 
have A(x)w = ReA^ = sin(#)/(l + cos0). For an ide- 
ally strong measurement, instead, (W — > 00) A(x)s — 
2ReA w /(l + \A W \ 2 ) = sin(6). For the allegedly optimal 
preparation, A(x)o = 1/A(x)s = 1/ sin(#). In the range 
tt/2 < < it, the inequality A{x)s < A(x)o < A(x)w 
holds, showing that the probe preparation is not optimal. 

We now proceed to show how to obtain an arbitrary 
shift by properly preparing the probe. We consider the 
functions 
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with support in [— mr/2, rwr/2], n G N, where a is an 
arbitrary real number and G(k) is a primitive of |-B(A;)| -2 , 
i.e. G'(k) = \B(k)\~ 2 . For n — > 00 the function is not 
normalizable, but gives rise to a finite shift. The initial 
and final shift of the probe are 



and 
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so that 
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We notice that the coefficient of a is zero only in the 
trivial case A w = ±1 (when the shift is A(x) = ±1 for 
any initial preparation of the probe, see the analysis be- 
low for a discussion of this case), otherwise it is always 
strictly positive. Hence, by properly choosing a, A(x) 
can be made equal to any value. 

A careful analysis reveals the following errors in 
Ref. Q: 

1. As Ref. [l| is implicitly looking for an extremal of 
the action (which is called a Lagrangian therein) for 
fixed values of £i(fc) at fc_,fc + , and these extrema 
of integration may be finite, Eq. (SSH23) and the 
first line of Eq. (SSH29) are not in general the 
correct definitions of the averages. They should be 
replaced by the real expression 
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with £ = B^i and £ = £j, respectively. Conse- 
quently, the reasoning that fixes the interval of in- 
tegration to [— 7r/2,7r/2] does not hold. However, 
since |£i| 2 is periodic with period 7r, this error has 
no consequence for the result. 

2. The integrand contains ^ which has a divergence 
at k = ±7r/2 for the wave- function given in Eq. 
(SSH15). As the <5-function given by contributes 
only a real term to £(fc)*£(fc)', this error has no 
consequences as well. 

3. Reference [l[ finds an extremum for (i)/, not for 
(x) f — (x)i. True, one can shift the initial wave- 
function in order to have (x)i — 0, however the 
variation is being made over all wave-functions. As 
shown below, this error is the cause of the incorrect 
claim. 



Appendix A: Correct treatment of the variational 
problem 

The wave-function of the probe after the measurement 



is 



(k\f) = if{k)/N) 12 = B{k%{k)/N^ < A I ) 

with the normalization 
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We shall consider that £i(fc) may not be normalized, 
so that the initial normalized wave-function is (k\i) = 

1/2 

£ i i(k)/N i with the normalization 
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The quantity to be maximized is the pointer shift (from 
now on we omit the functional dependence on k) 



A(x) =(f\x\f) - (i\x\i) 
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The problem is well posed when either the integration is 
over a finite interval fc_ , fc + with £j fixed at the extrema, 
or when the integration is over the real axis and it is 
required that limfe_ > .± 00 £j = 0, with the vanishing of £j 
fast enough to guarantee that N is finite. 

Since the premises of Ref. [J appear problematic, let 
us rederive the extremal conditions for A(i). We start 
by noticing that the action A(i) is invariant under the 
operation £i(fc) — > Ce ix ° k ^i(k) for arbitrary C £ C, xq E 
M. Thus the solution to the extremal problem, if it exists, 
is not unique. Physically, this symmetry means simply 
that shifting the initial wave-function £,i(x) by Xq does 
not change the shift of the pointer. By varying the action, 
we get 
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Hence 
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with \B\ 2 = Nf/Ni. We notice that Eq. jA6|) differs from 
Eq. (SSH25), since the extra term originating from — (x)i 
in the action is missing in Ref. 

A special case is A w — ±1, which gives B = exp 
Then the coefficient of ^ in Eq. (|A6j) vanishes identi- 
cally, and the acceptable solution is A(x) = ±1, with £j 
arbitrary. This has a simple physical interpretation: if 
B(k) — exp[=Fifc], the effect of the interaction (which, 
we remind the reader, is of multiplying the initial wave- 
function by B) in the position representation is to shift 
the wave- function by ±1, so that for any initial state, the 
shift is always the same. 

For A w 7^ ±1, we have a first order equation, with 
singular points (the continuous function \B\ 2 certainly 
equals its average for some k). The quantities appearing 
in the equation are 
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and 



We have 



B*B' = — - (1- \A W \ 2 ) sin(2fc) + Im(A„) cos (2fc) 



iRe(A w ). 
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rB n = d\B(k)\ 2 

equation Eq. (|A6[) is hence rewritten as 



We remark that 2Re(g* B') = a| ^' The differential 
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with 13 = LB r — \B\ 2 . Its solution is readily found to be 



exp < —i(x) ff 



A(x)|B| 2 -Re(A^)] ff(fc)} 
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with -ff(fc) a primitive of -D(fc) 1 ■ 

The parameter |i?| 2 must be found self-consistently by 
plugging Eq. (IA10p into 
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Thus the integration region [£;_, must include the ze- 
ros of D(k), and this implies that £i(k) is not normaliz- 
able. 

The average A(x), however, is well defined even if £j is 
not normalizable. By plugging Eq. (|A10[) into Eq. (|A4j) 
we have an identity. Thus, all the functions given in 
Eq. (|A10[) provide an arbitrary value A(x) for the shift. 

In the main text, see Eq. ([3]), we have provided a dif- 
ferent family of normalizable functions that provide an 
arbitrary A(x) as well. 
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